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ABSTRACT 

Hydrogen is injected from a downstream facing step in a wall into a high- 
temperature stream. Temperature and hydroxyl radical concentration are meas- 
ured downstream of the injection plane by means of ultraviolet absorption spec- 
troscopy. The experimental results are compared^ with theory which is based on 
a finite-difference solution of boundary-layer equations. Finite-rate kinetics 
equations are included in the analysis. The analytic predictions are also com- 
pared with previously obtained experimental results which are based on probe 
measurements. Comparison is made between calculated and observed ignition dis- 
tances . 


NOMENCLATURE 

A,B,C,D different chemical species 

A pre-exponential constant in the rate expression 

[A] mole fraction of A 

Cp specific heat 

E activation energy 

t reaction rate 

H total enthalpy 

h static enthalpy 

k specie number 

A reaction number 

M third body 

"''Associate member. 
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number of species 

number of reactions 

coordinate normal to the streamline 

temperature exponents in the rate expression 

static pressure 

magnitude of velocity 

gas constant 

universal gas constant 

streamline coordinate 

static temperature 

x component of velocity 

any dependent variable 

molecular weight 

Cartesian coordinates 

molar concentration 

rate of change of 

third body efficiency 

mass fraction 


characteristic length in the viscosity model 
flow angle 

effective turbulent viscosity 
stoichiometry coefficients 
density 

stream function 
net specie exchange rate 
Subscripts : 

A , B ,C,D different chemical species 
a,b respectively, forward and reverse reaction 

i,j grid points (Fig. 8) 

k specie number 

i reaction number 
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M third body 

<° main stream 


INTRODUCTION 

Recently, analytical modeling of the diffusive combustion and ignition 
processes (1-3) have received attention due to their application to the design 
of the supersonic-combustion ram- jet engine. 

While most of previous investigations were concerned with axi symmetric 
configuration, this paper considers a two-dimensional configuration in which 
hydrogen is injected from a rearward facing step in the direction parallel to 
the main stream. 

In the literature dealing with similar reactive-flow problems, it is 
usually assumed in the analysis that the chemical source terms may be modeled 
by expressing the reaction rates in terms of time-average quantities. For ex- 
ample, this assumption was made in reference (1) where individual combustion 
reactions were considered, and in references (2_,3) where global expressions 
were used to model chemical reaction. It is also usually assumed that there is 
no direct influence of chemical reaction on the turbulent transport mechanism, 
i.e. , that the turbulent transport model derived from a similar but nonreactive 
flow field applies. Some degree of success In predicting experimental data has 
been achieved (1-3) by using these assumptions. 

The present study provides an additional test for these assumptions by 
comparing the experimental results with predictions. 

The existing experimental information on reactive flows is mostly related 
to the axisymmetric free-jet geometry so that presented experimental results 
extend the range of available experimental data. 

The experimental conditions in the main stream corresponded to Mach 
number 2.44, static temperature 1270 K, and static pressure approximately equal 
to one atmosphere. Conditions in the hydrogen jet corresponded to sonic ve- 
locity, room total temperature and static pressure equal to that in the main 
stream. 


EXPERIMENT 
Test Section 

The test section which was used in the experiment is the same as in refer- 
ence (4) ; it is illustrated in Fig. 1. For the OH absorption measurement, the 
downstream windows were replaced with brass plates containing a line of 0,32 cm 
diameter holes, thus eliminating the windows as a source of beam attenuation. 
Two sets of plates were fabricated, enabling measurements to be made on 0.32 cm 
centers across the stream, 22.8 cm downstream of the injection step. A more 
detailed description of the test hardware Is presented in reference (4) , 

Optical System 

The optical technique used to determine OH concentration and temperature 
In this study has been developed and used in references (5^7) . This method 
relies on absorption of the narrow lines emitted from a hydroxyl radical source 
by broader lines in the combustion gas. 

Figure 2 illustrates the arrangement of the optical system. The OH line 
spectrum source was a water-cooled , end-view capillary discharge tube (6^2_) * 

The emission from the source was chopped and focused on the test section cen- 
terline with a quartz lens. The beam emerging from the other side of the test 
section was reflected by a spherical mirror to a focus on the entrance slit of 
1/2 meter grating monochrometer which was used in the second order, thus pro- 
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viding a good separation of the spectral lines of interest* The spectral lines 
were in the 0-0 band starting around 0*307 pm, The line intensity was measured 
by a photomultiplier located at the exit slit. The photomultiplier output was 
amplified with a lock-in amplifier and the signal recorded on an x-y recorder. 

Procedure 


The high- temperature gas stream was produced by burning a hydrogen- 
nitrogen gas mixture with liquid oxygen at high pressure. The oxygen content 
of the products was about 21% by volume* The Mach number in the gas stream was 
2*44, the static temperature 1270 K, and static pressure about atmospheric* 
Hydrogen was injected at sonic velocity, room total temperature, and approxi- 
mately atmospheric static pressure. 

Because of the use of heat absorbing hardware, the run times were limited 
to about 3 sec. In order to determine the temperature, absorption measurements 
for a number of different spectral lines are required. It was /necessary, 
therefore, to duplicate the same run conditions a large number of times, chang- 
ing the monochrometer setting between runs. A reading of emission line inten- 
sity with no absorbing gas present was recorded before and after each run. If 
the pre- and post-run levels were not the same, within about 2 percent, the 
data from the particular run were discarded. For consistency, the recording 
trace was always read at the same point relative to the end of the run* 


After a series of runs was completed at one location, the optical system 
was realigned through the next hole and another series of runs initiated. 

Overall, measurements were taken at 9 different locations across the 
stream, ranging from 0,6 to 6.0 cm from the wall, all at the axial location 
22.8 cm downstream of the hydrogen injection step* 

A brief description of the data analysis procedure is given in Appendix A. 

ANALYSIS 


The Equations 

It is customary in the numerical solution of jet mixing problems to re- 
place the normal coordinate y with the stream function by using the fol- 
lowing transformation. 


dy 


pu 


( 1 ) 


For the wall-jet problem this system becomes less suitable, because at the 
wall, the flow rate approaches zero and for the constant Increment AiJj, Ay 
becomes very large. This is just the opposite of the desired variation of Ay. 
However, in the streamline coordinate system (s,n) the boundary-layer equa- 
tions are as simple as in the (x,^) system and at the wall An *» Ay. It is 
noted that in the numerical procedure An must be allowed to vary, which ne- 
cessitates the formulation of the problem for a nonuniform grid* However, 
this does not complicate the solution since in the numerical calculation of 
the boundary -layer , the use of nonuniform grid is necessary for obtaining ac- 
curacy and convergence of the solution (8-10) * 

Assuming that turbulent Prandtl and Lewis number are unity, the time aver- 
age boundary layer equations written in terms of effective turbulent viscosity 
P are. 
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and the equation of state is 
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A \ 

p = pTR o J WT 


(5) 


( 6 ) 


The assumption that Prandtl and Lewis numbers are equal to one has been 
made because the equations are much simpler in this case and it appears that 
there is no concensus in the literature as to their values. 


The Effective Viscosity 

The turbulent viscosity formulation is as in reference (;4) . It is essen- 
tially based on the Herring -Mel lor ( 11) kinematic eddy viscosity model devel- 
oped for a turbulent boundary layer* The modification of the model is in the 
near region where the free-shear layer has not yet merged with the developing 
wall boundary-layer. In this region, the kinematic viscosity in the free- 
shear layer - analogous to the viscosity in the outer^portion of the boundary 
layer ( 11 ) - Is taken to be proportional to length 6^ given by 
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where y 3 and correspond to free-shear layer boundaries, Fig. 3, deter- 

mined so that the deviation from uniform flow on either side is about 0.1 per- 
cent. All constants in the model are as determined originally in refer- 
ence (11) ♦ 

Using this model in reference (4) good agreement between experimental and 
computed composition profiles was obtained for the case of no reaction. It is 
noted, however, that the model has not been tested by applying it to other ex- 
perimental data. 


Chemical Species Production Terms 

All the reactions considered in the analysis could either be represented 
by a bimolecular reaction 

v a U)A + v B (fc)B t v' c WC + Vp(R)D (7) 

or a three-body recombination reaction 

v a U)A + M t Vg(Z)B + v^U)C + M (8) 

where i = A,B,C and D are the stoichiometric coefficients for reac- 

tion number £, and M is the third body. Then, for example, for specie A, 
the chemical production rate is given by 

NR 

*a = Y2i i ^ w ■ < 9 > 
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where 


v^(&) \>b(£) 

y b 


Vq(1) v^U) 


f b^ Y C Y D 


for reaction (7) and 

v A (Jt) vi(S,) Vr(£) 

“l - f a U)Y A Y M W ' f b WY B Y C V*> 

for reaction (8)* The molar concentration Y M (Jt) depends on the assigned third 
body efficiencies k * 1* * • *> N 


V*> V 


k Y k 


( 10 ) 


k^l 


The temperature dependence for the forward reaction rate f a , for example, is 
given by 


— n a (Jt) 

f U) = A (A)T exp[-E U)/RT] (11) 

a a a 

Eighteen chemical reactions and nine chemical species were considered in the 
analysis. The reaction scheme and the rate information are essentially the 
same as presented in reference (12). 

Numerical Methods 

The system of equations (2) to (6) , together with appropriate equations 
for the production of species is a quasi-linear system of parabolic equations. 
Since the chemical specie production terms do not involve derivatives, they 
could be computed from values of dependent variables calculated from the pre- 
vious integration step. However, because of the well known difficulty in in- 
tegration of the one-dimensional flow problems involving chemical kinetics 
equations, special techniques were used in the past to evaluate the effect of 
chemical production terms. 

In reference (1) diffusion and reaction steps were solved independently 
and then the respective solutions were superimposed at each integration step. 

A special integration technique which had been developed for the one- 
dimensional problem was applied for the reaction step. This technique is 
based on linearization of the chemical production terms and the solution of 
the resulting set of linear differential equations by the use of parabolic 
fits for the specie concentrations. For the solution of the diffusive step, 
the explicit difference scheme was applied. 

Because of the presence of the boundary layer in the wall -jet problem 
considered in this report, the integration of the diffusion step by the ex- 
plicit method becomes impractical and so an implicit integration method is 
used. For the integration of the reaction step, the method of Lomax and 
Bailey (13) has been selected. 

These two steps and the coupling scheme are discussed more fully in 
Appendix B . 


The Initial Conditions 

The initial composition in the main stream was calculated by using the 
general one -dimensional computer program of Bittker and Scullin (12_) . The 
calculation was started in the gas generator and terminated at the test sec- 
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tion entrance. In addition to the chemical reactions involving hydrogen- 
oxygen species , several reactions involving the NO radical were also included 
in the calculation. However, the inclusion of these reactions had no effect 
on the ignition delay on the one-dimensional basis (4) and, therefore, they 
were omitted in subsequent two-dimensional calculations. 

The calculated initial composition in mole-fractions is [H 2 ] " 0.193x10“^, 
[0 2 ] = 0.208253, [H 2 0] = 0.3410570, [H] = 0.8xl0 -6 , [ 0 ] =■ 0.29x10-4, 

[OH] - 0.3419x10-3, [H0 2 ] = 0.48xl0" 5 , [H 2 0 2 ] = 0.27xl(T 5 , and [N 2 ] = 0.4502913. 

The calculated static temperature in the main stream proved to be somewhat 
higher than the experimental value. This discrepancy is believed to be due to 
the heat transfer effects which have not been included in the computation. 

Since it was essential to reproduce the static temperature correctly, the com- 
position of the reactants in the gas generator chamber was adjusted by reducing 
the hydrogen content slightly until the resulting static temperature was suffi- 
ciently close to the experimental value. 

The initial velocity and temperature distribution determined on the basis 
of measurements are as presented in reference (4). 

In the experiment, largely because of the boundary-layer build-up on the 
test section walls, a small change in the mainstream flow variables was noticed 
between test section entrance and exit. In the analysis, this change was as- 
sumed to take place linearly. The equations for the main stream velocity, 
static temperature, and static pressure are: 

- 1791 - (x/ 35 .6)128. 8 m/s (12) 

T^ - 1259 + (x/35. 6)76.9 K (13) 

= 9.24xl0 4 + (x/35.6)2.48xl0 4 N/m 2 (14) 

COMPARISON OF ANALYSIS AND EXPERIMENT 
Ignition Distance 

The measured ignition distance determined from the boundary of the lumi- 
nous zone was on the average about 15 cm, whereas the calculated ignition dis- 
tance determined on the basis of a 5 percent rise in static temperature above 
the mainstream value was 16.5 cm. 

The agreement between measured and calculated values is satisfactory, 
particularly in view of the large temperature difference between the hydrogen 
jet and the main stream. 

Comparison of Predicted Profiles with Optical Measurements 

The static temperature profiles at x = 22,8 cm are presented in Fig. 4, 
while the OH mole fractions are presented in Fig. 5. The calculated profiles 
are presented for two sets of main-stream composition. One set corresponds to 
the calculated composition. In the other, the OH mole fraction is adjusted to 
correspond to the measured main-stream value. The mole fraction of other rad- 
icals were also adjusted in the same proportion. 

Good agreement is noted in Fig. 4 between predicted and experimental 
values of maximum temperature, when the prediction is based on experimental 
value of OH in the main stream. However, predicted temperature profile is 
considerably narrower. A much more gradual approach to the main-stream tem- 
perature is noted in the experimental profile. 

For the case of calculated main-stream composition the predicted maximum 
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temperature is about 200 K higher than the experimental* However, the average 
width of the calculated temperature profile now more closely corresponds to the 
experimental value. 

The OH profiles in Fig* 5 indicate considerable differences between theo- 
retical and experimental values. The magnitude of the deviation appears to be 
proportional to the OH level which is the opposite of the trend observed in 
Fig. 4 for temperature profiles. Qualitatively, this could be explained by the 
way these two measurements are influenced by the nonun if orraities across the 
optical path length introduced by turbulence and by the presence of the test- 
section side boundary-layers. The hydroxyl radical distribution is averaged 
evenly, independently of the magnitude, whereas, in the case of temperature, 
the indicated value more closely corresponds to the maximum across the path 
length. 

The optical path length used in the data reduction has not been corrected 
for the presence of the side boundary-layers. Because of associated effects 
due to test-section corners, the correction is probably more significant in the 
region close to the wall (small y) . In the main stream, reasonable agreement 
is noted between experimental and predicted OH mole fractions. 

To illustrate to what extent the hydroxyl radical concentration varies in 
the x direction in this region, also plotted in Fig. 5 is the calculated pro- 
file of x = 20. 1 cm. It is seen that in a distance of less than 3 cm the OH 
mole fraction varies by more than an order of magnitude. 

Comparison of Predicted Profiles with Probe Measurements of Ref. (4) 

In reference (4) the composition and the flow fields derived from probe 
measurements at x - 35.6 era were compared with corresponding predicted pro- 
files which were computed by assuming that combustion reactions are in chemi- 
cal equilibrium. The composition profiles from reference (4) are reproduced 
in Fig. 6. It can be seen from Fig. 7 that the agreement between theory and 
experiment is now better and that, therefore, as anticipated, some of the ob-^ 
served shift of the calculated profiles in relation to the experimental in 
Fig. 6 is mainly due to the equilibrium assumption* 

In Fig* 7, the most noticeable departures of computed profiles in rela- 
tion to experimental profiles are: (1) lower mole fractions of hydrogen and 

higher mole fractions of water vapor in the region where y < 2 cm, and 
(2) more gradual profile slopes in the region where y > 2 cm. The latter In- 
dicates that the assumed mixing rate in the outer portion of the boundary layer 
is too high. It is also possible that a better prediction of the experimental 
results could be obtained in the region where y < 2 cm by assuming a lower 
mixing rate in this region, because this assumption would cause an increase in 
hydrogen concentration and a reduction In water-vapor concentration. This is 
evidenced by comparing Figs. 6 and 7, The thickness of the mixing /boundary- 
layer region in the predicted profile in Fig. 6 is smaller than in Fig. 7. 
Because viscosity is proportional to this thickness, it follows that viscosity 
is also smaller , and it can be seen that the hydrogen mole fraction at the wall 
is predicted better in Fig. 6. 

Taking into account the results of the previous section and of the fact 
that the ignition distance Is fairly well predicted, the overall effect con- 
cerning the rate of reaction is that it is fairly well modeled up to the igni- 
tion point and that it is somewhat overpredicted downstream of it. 

Comparison was also made between predicted and experimental profiles of 
total temperature and the Mach number. The total temperatures agreed fairly 
well and markedly better than in reference (40 . However, significant differ- 
ences, particularly in the region close to the wall were noted in the case of 
Mach numbers. The difference, it is believed, arises because of excessive 
flow deceleration in the test section caused by boundary-layer effects. A 
similar conclusion was made in reference (4) comparing Mach numbers for the 
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case of no reaction. 


CONCLUSION 

For the wall-jet configuration, the use of the time-average concentration 
and the molecular reaction rate expressions together with a kinematic eddy- 
viscosity representation of the turbulent transport phenomena results in (1) a 
reasonably accurate prediction of the ignition distance, (2) a fair prediction 
of the magnitude and the relative position of the peak static temperature in 
the early stages of the flame development, provided that the hydroxyl radical 
level in the main stream is based on the experimental value, and (3) a fair 
prediction of the stable specie distribution in the diffusion-reaction zone. 

The departure of theoretical predictions from experimental results indi- 
cates that the assumed turbulent transport model causes the rate of mixing in 
the outer portion of the mixing-boundary layer to be too high and that the re- 
sulting theoretical rate of reaction is somewhat over-estimated downstream of 
the ignition point. 
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APPENDIX A 

The absorption coefficient, a, is given in terms of the path length, L, 
the measured line intensity with no absorber present, I Q , and the line inten- 
sity after passing through the absorbing gas, I, 

a ° ln(l 0 /I)/L (15) 

On assuming rotational equilibrium, temperature is determined from the 
Boltzmann plot, an example of which is illustrated in Fig. 8. The ordinate in 
this figure is the natural logarithm of the ratio of absorption coefficient to 
the product of the relative transition probability, A^, and a correction fac- 
tor for vibrational-rotational interaction, Tj'j" (14) . The abscissa is the 
wave number of the rotational level, divided by the Boltzmann constant k 
and multiplied by the product of Planck's constant h, and the speed of light, 
c. The slope of the line in Fig. 8 is the inverse of the temperature. 


The number density of the OH radicals, Nqh, is given as: 


N 


OH 


2.4x10 


12 Q R Q V b D° 
FA k T J'J 


II 



(16) 


where and Qy are, respectively, rotational and vibrational partition 

functions; is the Doppler half-width of the observed spectral line; and F 
is a constant taken from reference (6). 

APPENDIX B 


Integration of the Diffusion Step 


In the absence of reaction the problem reduces to the solution of the 
turbulent boundary-layer coupled with the jet-mixing flow. The differenced 
form of the boundary-layer equations (2) to (4) for the mesh defined in Fig. -9 
is given by 


(pq) i,j 


V - v 


As 


An l,j-l/2 \ ^ +1/2 


V - V 

i+l,j+l i+l,J 

An . 


- Vi 




V - V 

i+l>J 1+1, .1-1 

An. , 
j-i 


(17) 


where V stands for either q, H, or a^.. 
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The continuity equation (5) is integrated by centering 0 halfway between 
the grid points. 


(pq) 


i+i , j - (pq) i, j + 


As . 
1 




6 i+l ,.1+1/2 


- 0 


An 


i+l.j-l/2 _ Q 


(18) 




Once 6 is known new grid spacings can be computed by using simple geometrical 
relations. For example, Ani+2. is given by 


An i+l,j-l/2 " in i,j-l/2 + As i,j (e i, j+1/2 ‘ 0 i,j-l/2 ) (19) 

Because the difference equations are implicit* it is noted that the nonuniform 
grid spacing does not introduce any stability problems. 

Because of large variation in viscosity, it was found necessary to vary 
the initial An spacing. This was done so that approximately 

2 

( P q / U ) An - constant 

without allowing an excessively large rate of increase of An away from the 
wall. In practice* following this guideline usually resulted in a mesh with 
several constant An spacings next to the wall followed by a section where An 
increased by a factor of 1.5 to 2 after which An was again constant. Total 
number of points in the normal direction was typically about 30. 

For the jet-mixing problem, the results obtained using this integration 
scheme were in agreement with results obtained by using a program which em- 
ployed the stream function transformation (1) and a uniform grid spacing. 

In the case of a turbulent boundary layer, the present program predicted 
somewhat lower values for the displacement thickness in comparison with the 
calculation of Herring and Mellor (11) . The normalized velocity profiles, 
however, agreed very well. In this regard, it may be noted that slight dis- 
crepancies in the prediction of turbulent boundary layers which apparently re- 
sult from the use of different numerical schemes have been observed in the 
past (15) . 


The Reaction Step and the Coupling Scheme 

Since the effects of the chemical reaction are evaluated independently of 
the diffusion process, equations (4) and (3) reduce to 


dy, 

(l/pq)Y k (T,p,Y l Yjj) i k = 1, • • -,N (20) 



W C p> 




( 21 ) 


where 


Y k - V w k 

The system of equations is completed with 
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( 22 ) 


1 11 
which is obtained from equation (6) by differentiation. 

The integration method is based on an implicit difference scheme (13) 
which is applied to the linearized system of equations. Linearization is car- 
ried out uniformly with respect to all the dependent variables 
T, and p. The derivatives are replaced using the simple formula 

* &)J 

where V represents any dependent variable. Substitution of (dV/ds)-[+i and 
(dV/ds)j_ in (23) is performed by using equations (20) to (22) upon lineariza- 
tion. It is noted that information at only one point (i) is needed in order to 
compute variables for a new point (i + 1) . 

The computation scheme differs from that in reference ( 13 ) in that the co- 
efficients in the linearized equations are computed from derived expressions 
rather than by employing numerical differentiation. This improves the effi- 
ciency of computation since the number of times such functions as those de- 
scribing the rate of reactions are evaluated can be significantly reduced. 

The linearization has been carried out for a general case of arbitrary number 
of reactions. 

On the one-dimensional basis, the results obtained by this method were in 
good agreement with Morretti f s solution (16) which is based on the exact solu- 
tion of the linearized system of differential equations. 

The coupling of diffusion and reaction steps can be accomplished in sev- 
eral ways. In the present program it is carried out by first solving the dif- 
fusion equations and then the reaction equations. The values of Yj_ , . . . ,Yfl, 
p, and T calculated at the end of diffusion step are used as initial values 
for the reaction step. The spacing As is controlled by the reaction step. 
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Figure 2. - Schematic of optical system used to measure temperature and OH 
concentration by UV absorption. 
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Figure 3. - Schematic of the near flow-field. 
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Figure 4. -Temperature profiles at x = 22.8cm. 
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Figures. - Hydroxyl radical profiles. 
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Figure 6. -Theoretical and experimental composition profiles. 
Equilibrium chemistry. (From reference [4]) 



Figure 7. - Theoretical and experimental composition profiles. 
Nonequilibrium chemistry. 
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Figure 8. - Boltzmann plot of OH line absorption measure- 
ments 1. 85 cm from wall. Symbols in the plot denote 
particular rotational lines. 
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